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Abstract 

Pair correlations and the magnetic susceptibility of electrons in a spherical cavity are 
studied both for grand canonical and the canonical ensemble. The coupling constant of 
the BCS Hamiltonian is adjusted to experimental values of the gap parameter. The gap 
parameter is found to increase for small grains as a consequence of the pronounced shell 
structure in the spectrum of the spherical cavity. The sharp phase transition at is 
smeared out for the canonical ensemble. The strong paramagnetic susceptibility of the 
normal electrons in the cavity is reduced by the superconductivity, but it remains positive. 
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1 Introduction 



Whereas many properties of superconducting bulk metals at low temperatures have been 
well studied, the influence of size effects on the same properties of small metal grains or 
clusters is still a problem of a current interest. The minimal size of grains in which the 
phenomenon of superconductivity can be observed has not been not established by now. It 
has been known for about three decades that metal grains on the nanometer-scale display 
superconductivity 

The focus of the paper are the electronic finite size effects in the grains, which reveal 
themselves not only in the discreteness of the energies of the electron states but also in 
their shell structure. This means, that the level density strongly fluctuates as function of 
the energy. In regions with high level density the pair correlations are enhanced. This 
kind of enhancement has been demonstrated for electrons confined to a slab and to a 
cube P] . In this paper we are going to investigate the same effect as well as the magnetic 
response of small grains by means of yet another very simple model. It is assumed that 
the electrons are confined inside a sphere or hemisphere with a perfect surface and that 
only the electrons on the highly degenerated Fermi level take part in the pair correlations. 
The spherical model has a very pronounced shell structure. In this respect it contrasts the 
extended Landau - Ginzburg model for grains [Q, which does not take into account any 
shell structure. Comparing these two extremes will provide insight into the properties of 
real grains with electronic shell structure. Though the model is a very strong idealization 
of realistic grains it has the advantage that the exact solutions of the many body problem 
with finite particle number are known. This permits to exactly calculate finite size effects 
due the conservation of the electron number as well as the enhancement of the magnetic 
response in a finite system. It also allows to judge the typical approximations (mean field 
approximation, grand canonical ensemble) one has to resort to in more realistic models of 
the electron system of the grain. In particular, it will be demonstrated that the strong shell 
structure of the (hemi) sphere generates a dramatic enhancement of the pair correlations 
with decreasing size, resulting in a growth of the transition temperature. 

The experimental transition temperatures Tc{N) for grains with N delocalized electrons 
turn out to be larger than or equal to the bulk value. For Al, a weak-coupling supercon- 
ducting metal, the measured ratio / = Tc{N)/Tc{bulk) grows considerably with decreasing 
N: f ^ 1.5 for ~ 10^ and / ^ 3 for ~ 10'^ S 0- Grains of intermediate-coupling 
metals show a more moderate growth of Tc{N): f ^ 1.2 for In grains with A^ ~ 10^ |Q, 
/ 1.1 for Sn grains with A^ ~ 10^ For strong-coupling Pb the transition tempera- 
ture Tc{N) is about equal to the bulk value down to grains with A^ ~ 10^ 0]. Recently, 
Black, Ralf and Tinkham |]^, [1^ provided additional evidence. By tunneling experiments 



on grains of i? ~ (5-10) nm, the pairing gap A has been measured to be as large as ~ 2 
times the bulk value. 

The increase of the transition temperature and the pairing gap has been attributed to 
an increase of the coupling constant A with decreasing size of the grains. Various theoretical 
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models have been suggested for explanation of such an increase [|], ^ |1^, |T^, |T^, |T^, |T6|, |17 . 
In this work we do not investigate the possible mechanisms behind such a change of the 
coupling constant. Rather, we choose a phenomenological approach, introducing an explicit 
dependence of the effective interaction strength A on the electron number A^. We are going 
to demonstrate that in studying the N - dependence of the coupling constant the shell 
structure has to be taken into account. 

The BCS theory of our model is exposed in sections ^ and ^ We limit ourselves to Al 
grains with 10^ < N < 10^, for which the best experimental data on A are available. In 
section |., these are used to determine G{N) by comparing the empirical and calculated 
values of Tc{N). 

It is known |]|, |T8| that the specific heat C(T) of superconducting metallic grains does 
not show the singularity at seen in bulk metals. For grains, the function C{T) has a 
maximum at Tq that is smaller than Tc for bulk metals. With decreasing grain size this 
maximum shifts to smaller temperature and the width of the peak becomes wider. The 
question arises, how to define a transition temperature Tc from the non-singular function 
C(T). Calculations in the grand canonical BCS approximation cannot reproduce the ob- 
served T-dependence of C. Therefore, in section |^ the infiuence of the superconductivity 
on C(T) is studied within the canonical ensemble. The results suggest a definition of Tc, 
which coincides with the one in BCS theory as well as with previous definitions, introduced 
for interpreting the measurements of C{T) p, |T8| and electromagnetic properties [T8[ 
in grains. 

Superconducting bulk metals show the Meissner effect. The diamagnetic susceptibility 
takes its maximum, compensating completely external magnetic field inside the metal. For 
grains with sizes comparable or smaller than the coherence and penetration lengths an 
incomplete compensation is expected, which will be investigated in section ^. The issue 
is complicated by the fact, that in normal (non-superconducting) grains the susceptibility 
at low temperatures strongly deviates from its bulk values. Both para- and diamagnetic 
enhancements appear a certain electron numbers, refiecting the electronic shell structure 



of the grain. This has recently been studied in M, EOl, ET] , where the references to further 



work can be found. Hence, an intricate interplay between the paramagnetism due to 
the shell structure and the diamagnetism due to the superconductivity is expected for 
superconducting grains, which will also be studied in section |^. 



2 The spherical cavity model in BCS approximation 

We consider electrons in a cavity of radius R = tqN^^^ , i. e. we employ the spherically 
symmetric rectangular well with infinite walls as a simplified model of the mean field that 
confines the delocalized electrons of the grain. The energies et are given by the roots of 
the Bessel functions. 

Thus, we describe the system of A^ electrons by means of the effective electron Hamil- 
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tonian 

H - fiN = ^(ej - fi)afat - G{N) ^ afafasas. (1) 

The interaction strength G is treated as a parameter that is fixed by comparison with the 
experimental pair gaps A(A^) measured for different grain sizes. The pairing interaction 
in (1) acts only among electron levels inside the interval \et — ^\ < Tiujd-, where the Debye 
frequency is taken to be the same as in the bulk. 

The treatment of the Hamiltonian in BCS mean field approximation is a standard 
problem in nuclear physics. For example it is exposed in the textbook [^]. We have 



solved the problem numerically. If a grain contains 10^ or more delocalized electrons the 
single-particle spectrum posses a high degeneracy in the vicinity of the Fermi level [F). 
The average orbital momentum / is of order N^^^. Each level can be occupied by about 
4 ■ N"^/^ electrons. The number of these degenerated levels inside the Debye interval 
varies from 1 (A^ ~ 10^) to 5 (A^ ~ 10^). However, as demonstrated in fig. 1, practically 
it is sufficient to take into account the Fermi level alone because for particle numbers 
(10^ < N < 10^) the distance to next levels is always much larger than the pairing gap 
A, found by means of the full BCS equations. The isolation the Fermi level and its 



high degeneracy are the conditions to apply the single shell model |Q for pairing. This 
approximation takes into account only the pair interaction among the electrons in the 
partially filled Fermi level. In the grand canonical BCS approximation, the single shell 



model was used in for a half filled shell. We consider the case of arbitrary numbers of 
particles in the shell and, in sections ^ and |^, the canonical ensemble of the exact many 
body states. 

In the single shell approximation, which we are going to study, the spherical cavity is 
equivalent with a hemisphere. This is seen as follows. The wave functions in the spherical 
cavity have good orbital angular momentum /. All wave functions with odd / are equal 
to zero in the the plane z = 0. Thus, they fulfill also the boundary condition for the 
hemisphere and are ( with the appropriate choice of the normalization) the wave function 
of this type of cavity. None of the conclusions drawn below depends on / being odd or even. 
Thus, they are valid for the hemisphere as well. Clusters on surfaces, the superconductivity 
of which is studied experimentally, have the shape of somewhat fiattened hemispheres. 

The single shell approximation leads to considerable simplifications because there are no 
summations over single particle states in the BCS equations and many quantities become 
analytical expressions of the temperature and the filling parameter 

A^ h 

n = —, M = 2/ + l, (2) 
2M' ' ^ ' 

where N^h is the actual and 2(2/ + 1) the maximal number of electrons on the Fermi level, 
which has M = 2/ + 1 magnetic substates and 2 orientations of the spin. There is only one 
quasiparticle energy 



E= V(e-/i)2 + A2, (3) 
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which is independent of n at zero temperature, 

m = (4) 

At finite temperature the quasiparticle energy is obtained as the sohition of the imphcit 
equation 

E = E(0)tanh{^}. (5) 

Independent of temperature, the difference between the Fermi level and chemical potential 
is given by 

e- fx = E{0){l-2n) (6) 
The pairing gap A has a maximum in the middle of the shell {n = 1/2). It is given by 



A = ^E^-ie- (7) 

Its zero temperature value is 



A(0) = 2E(0)^n(l -n). (8) 

As displayed in Fig. 2, a finite temperature T decreases A but its behavior as a function 
of n is similar to that at T = 0. Our calculations of A with all single particle levels inside 
the Debye interval give results practically indistinguishable from those in Fig. 2. This 
demonstrates the applicability of the single shell approximation for the considered values 
of N and T. 

The transition temperature Tc is defined as the temperature where A disappears. It is 
obtained from Eq.(|), substituting the value of the expression (j^) at A = for E, i. e. 

' ^' ^ ^exp{(|e-/i|/Tj} + l 

Taking into account the expression (j^) for the chemical potential the transition temperature 
is given by 

Un)=2E{0)r\n) , (10) 

1 + X 

^{n) = 2x~^ \n , x=\l-2n\, < a; < 1. 

1 — X 

As shown in Fig. 3, is symmetric with respect to exchanging n ^ 1 — n. In the middle 
of the shell (n = 1/2, a; = 0) the limit x — * in Eq. ([To|) gives the known result ||22| 



^{n = 1/2) = 4, which coincides with the value obtained Parameter for a cubic cavity. 

In the BCS theory of bulk metals there is the universal relation between Tc and 2A(0) 
which is the energy of the lowest two-quasiparticle excitation, 

2A(0) 2E(0) 

^^ = ^^ = 3.52. (11) 

-'■ c -'- c 
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In the single shell model at T=0 the two-quasiparticle energy 2E{0) is different from 2A(0). 
Accordingly, two different ratios can be considered, 



an) (12) 



and 



2A(0)/T, = 2eH^n(l-n). 
Both depend on the shell filling parameter n. 

3 Averaging over the shell structure 

The physical quantities we are interested in show rapid variations as functions of the 
electron number N, which are manifestations of the pronounced shell structure. We will 
average some quantities f{n) over an interval of particle numbers N, 



{f{y))av = I /(x)exp{— 



where x = N^^^ and 7 = 0.3. The reason for averaging consists in following. Firstly, the 
number of electrons is not exactly known and there is a considerable uncertainty in the grain 
radius. In measurement on probes containing many grains, there will be an experimental A^- 
distribution. Secondly, the phases of the fast iV-oscillations are sensitive to small deviations 
from the ideal spherical symmetry of the adopted model. In real grains this symmetry is 
certainly broken by the roughness of the surface due to the discrete ionic background, by 
deviations of the shape from a sphere or hemisphere and by impurities. These imperfections 
will shift and wash out the oscillations. Averaging out the shell structure oscillations has 
also been advocated in ref. [0 studying the superconductivity of electrons in a cubic cavity 
as well as in the studies of the enhancement of paramagnetism in mesoscopic systems (cf. 



e. g. PI). 

The averaged values of 2E{0)/Tc are shown in Fig.4. They are larger than 4 because Tc 
is less than Tc{n = 1/2) for n ^ 1/2 but E{0) is n - independent. The averaged ratios of 
2A(0)/Tc are close to the bulk value 3.52 of a weakly coupled superconductor. Since the 
same ratio is also found for the cubic cavity [|l^ it may be of general nature. 



4 Effective pairing strength in the BCS approxima- 
tion 

The two quasiparticle energy 2E{0) has been measured in refs. |P, [T^ for different A^ at 



sufficiently low temperature, such that the zero temperature expressions can be applied. 
Within our single shell model, eq. (|| ) relates them directly to the coupling constant G. In 
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order to reproduce the experimental - dependence of 2E{0) we must assume that G{N) 
deviates from being oc as in the bulk. The restricted number of data points does 
not allow a quite definite determination of this function. We have adjusted two different 
phenomenological expressions, 

G = gN~'', g = 1.94 meV, a = 0.47, (13) 

and 

G = gN~^ exp [-aiV-^] , (14) 
c/ = 3.21e1/, a = 25, /3 = 0.26 

Since often transition temperatures are measured, we show in figs. 5a and 6a the 
values of Tc obtained by means of (|T^) and the expressions (|TBp and (|I3p, respectively. For 
comparison, the values of 2i?(0) measured in refs. |P, |l^ are converted into "experimental" 
Tc values by means of eq. (IH), using the averaged values of ^. 

As seen in Fig.5a, expression (|13D is quite reasonable in the region of ~ 10'^ -r- 10^. 
Varying a one can obtain as an increasing (a < 1/3) or decreasing (a > 1/3) function 
of A^. In the single shell model the particular role of a = 1/3 is connected with the 
proportionality of Tc to GM i.e. to GN^^^ (M averaged is ~ N^^^). 

The A^ - independent bulk coupling constant is given by A = Gpp, where pp o(. N is 
the density of states at the Fermi- surface. Eq. ([I^) is constructed as a product of the bulk 
coupling constant and an A^ - dependent factor that accounts for the finite size effects. 
Accordingly, g is estimated using the coupling constant of bulk Al, A = gppN'^ = 0.4 and 
the Fermi-gas density pf{AI). Both in eq.(|l^ and eq.(|l^) there are two fit parameters. 



Eq.(|T4D gives the correct asymptotic behavior of G aX N ^ oo. The particular choice of the 
A^-dependence corresponds to decreasing G at small A^, which is reflected by the decrease 
of Tc and A at small A^ in Fig. 6. No physical significance is attributed to this decrease, 
because other choices of the A^-dependence of the factor are possible. 

The function G{N) obtained by fitting the spherical model to the data strongly deviates 
from the bulk-law G oc A^~^. It also deviates from the function G{N) one would obtain 
within the frame work of models that disregard the shell structure, like the one of In 
oder to illustrate this statement we assume equidistant levels near the Fermi surface with 
the spacing d = pp^. This model is exposed e. g. in |2^. For A^ > lO'^, when A{0)/d ^ 1, 



it is possible to replace the sums over the single electron levels by integrations and the well 
known expression 

E{0) = A(0) = hcuD/ sinh(^) ^ 2ujn ^M'-^) (15) 

for the bulk is obtained. Expressing the pairing constant as the product G{N) = d\f{N), 
the factor f{N) must grow with falling A^ in order to reproduce the experimental observed 
increase of A with decreasing A^. This is at variance with our fits. In the case of eq. (0), 
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f(N) = O.QN^-^^ and in the case of eq. (0), f{N) = exp(-25A^-°-26), which both decrease 
with decreasing N. 

Hence, if the shell structure in small grains is not taken into account, the effective 
coupling constant must be increased as compared to the bulk. On the other hand, in our 
spherical model the strong bunching of the electron already increases the value A(iV) such 
that the coupling constant must be attenuated in order to account for the more modest 
increase seen in experiment. However, we do not consider our fits as evidence for a decrease 
of the coupling constant. Due to its high symmetry, the spherical model over-accentuates 
the shell structure. Imperfections of different nature will make the level bunching in real 
grains less pronounced. The spherical model and the model with equidistant levels may 
be considered as the two limiting cases of maximal and no shell structure, respectively. 
The real grains lie somewhere in between. It is quite reasonable to assume that their shell 
structure is weak enough such that f{N) still increases with decreasing A^. 

In the context of the experimentally observed increase of Tc in thin films several effects 
have been discussed that lead to an increase of the coupling constant as a consequence of 
the reduced dimensionality |^, [Tl|, |12|, [T^, |14|, [T^, |16|, |T^. In particular, it has been pointed 
out that the surface phonon modes in films and small grains are expected to amplify the 
pair interaction. Our study shows that the electronic shell structure must be taken into 
account if these effects are quantitatively related to the experimentally observed increase 
of Tc, because the films often have a granular structure. Even for homogeneous films 
the quantization of the electron motion perpendicular to the surface must be taken into 
account, because ref. has demonstrated that the pair gap A is considerably larger than 
its bulk value (|15D when the thickness becomes as small as few times the Fermi wave length 
(assuming G = Xpp^)- 



5 Critical temperature in finite systems 



So far we have treated the pairing in the frame of the grand canonical ensemble. An import 
aspect of the small finite systems consists in the fact that the number of electrons in the 
grain is fixed and the the canonical ensemble must be applied. In ref. this question has 
been investigated on the basis of the the generalized Landau - Ginzburg equations. This 
approach does not take into account the discreteness of the electron levels and their shell 
structure. The spherical model in single shell approximation is simple enough to carry 
through the canonical statistics including the shell structure. 



As shown by Kerman pj], the pairing Hamiltonian can be expressed as the Casimir 
operator of the quasi spin algebra 



The eigenvalues E^^^ are characterized by the quasi spin Qg = (M — s)/2 or seniority. 
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which equals the number of unpaired particles in the shell, 



El^) = -G 



QsiQs + 1) 



In this section we consider only the particles on the Fermi level, i. e. = Ngt- Each 
eigenvalue (except the ground state one, Qo = M/2) is degenerated with the multiplicity 

The single shell model gives the energy of two quasiparticle excitations (states with 
s = 2) equal to 

) - i^f) = G [QoiQo + 1) - Q2iQ2 + 1)] = GM, 

i.e. it is exactly the same value (^) as in BCS, where GM = 2E{0). Thus, using 
the exact solution does not change the way to determine of G from the empirical two 
quasiparticle energies at T ~ [§, |lOl- However, data on in small grains require the 
consideration of the temperature dependence of the superconductivity in this model, which 
will be discussed below. 

All thermodynamical quantities can be calculated from the canonical partition function 
Z, which is 

N 

5:4exp{-(Ef)-EfVn, 

s=so 

where Sq = (sq = 1) and s are even (odd) integers if is even (odd). 

The absence of a sharp transition from superconducting to normal phase in finite sys- 
tems can be demonstrated, for example, by studying the temperature dependence of the 
internal energy < H >: 

< H, >= ^E^^^f ^exp{-(i?f ) -i^fVT}. (16) 
It is frequently represented in the form analogous to the BCS expression: 

< H, >= - Eex- (17) 

At r = 0, Acan = G{N + 2 \ A"^ \ N), where | A^) is the ground state function. It 
becomes A in the BCS approximation. E^x is analogous to the exchange energy in BCS 
( Eex{BCS) = —Gj^i^i)"^)- There is an arbitrariness in the choice of this term in (|l^. In 
ref. p6| it is chosen such that for half filled shell (A^ = M) it is equal to —CM/ A, which 
which is the temperature independent BCS value in the single shell model. The choice 

Epx =< Hn >It^oo= „ f — = — G — — \ (18) 

P M^oo 2(2M-1) ^ ^ 
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seems more appropriate to the problem, because for T — > oo, where the pairing disappears, 
< H > given by ( p!6D becomes equal to (p^S]). The difference between (0) and Eex{BCS) = 
—GN'^/AM is negligible at large M and A^. Inserting Ef.x given by (|T^ ) into the expression 
(|l|) for < Hp>, one gets for T = 



Acan(T = 0) =A(0)W1 + 



2M- 1' 



(19) 



where A(0) is the BCS parameter given by (^). Thus, Acan is larger than A(0) at T = 
but their ratio is practically equal to 1 for large M. The quantity Acan exceeds A{BCS) at 
any temperature, especially after where A{BCS) = 0. As seen in Fig. 7, it is smoothly 
decreasing and there is no sharp transition to the normal state. The larger M, the closer 
Acan approaches A{BCS). 

The specific heat C of the electrons, 

1 9^ Id 
C = —T^{TlnZ} = <H>, 
N dT^^ ^ NOT 

becomes in the single shell model 



C 



N{Ty 



^)'exp{-(Ef ) -^rVn- < H, >^ 



This equation can be represented in a G-independent form if we introduce the BCS critical 
temperature Tc as the unit of the temperature and energy, 

£^ = Ef)/T„ t = T/Tc, 

where is determined by (^) for given M and A^. As a consequence of and Eg being 
proportional to G, one has 



C 



k 



^J2dsel exp( 



t ■ 



1 
Z2 



J^ds^s exp(- 



t ' 



}• 



(20) 



In the grand canonical BCS approximation C can be represented as follows. 



C 



1 

N 



df~G 



In the case N = M the BCS gap equation (^) takes the form {t = T/Tc 



tanh(- 



A 1, 



A(0) 'A(0)t'' 
In the limits T — > and T Tc one obtains for the gap and specific heat: 

A 



t < 1 : 



A(0) 



l-2exp{-2/t}, ^^^ = ^exp{-2/t}. 
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A r- C(BCS) 

It is interesting to note that the jump in the C{BCS) at t = 1 coincides with the value of 



the Gorter-Casimir two-fluid model 25 



Fig. 8 shows C calculated by means of (pOl) . The case of the half filled shell N = M 
is shown for several values of M. For comparison, the BCS result is added. Our model 



reproduces qualitatively the T dependence of C observed in small grains |18|: when N 
decreases, the peak in C shifts to smaller t, it becomes lower and its width increases. A 
similar qualitative behavior of C was also obtained in the framework of the generalized 
Landau - Ginzburg model and the two- level model system The similarity of the 
T dependence of C found in the three models, which are very different with respect to the 
shell structure of the electron levels, indicates that the behavior of the function C{T) shown 
in Fig. 8 is a general consequence of the small and fixed number of particles. Fig. 8 indicates 
that at t = 1 (T = T^) the canonical specific heat in finite but not very small (M > 10) 
systems approximately attains the half of its maximum value. This observation correlates 
with empirical definition of Tc as the temperature at which a measured quantity takes 0.4 
or 0.5 of its maximal value. This definition was used in measurements of both electronic 



specific heats p, [T^ and electromagnetic quantities |jT^ |T8[. The critical temperatures 
determined this way practically coincide. Thus, the consideration of the theoretical specific 
heat in the single shell model canonical approach shows that the critical temperature 
calculated in the grand canonical BCS approximation can be adopted as the transition 
temperature for small systems. Our determination of the coupling constant G performed 
in the previous section on the basis of measured energies of two quasiparticle excitations 
and critical temperatures practically require no corrections. 



6 Magnetic susceptibility 

The magnetic susceptibility x of the electrons in the low field limit is given by the change 
of their thermodynamic potential 

Q{T, B) = {H-f,N + ujl{L. + 25,) + ^(x^ + y')) - TS (21) 
with respect to the external magnetic field B (assumed to be in direction of the z - axis), 

_ 1 d^n 

^ - ~VdB^^''='- ^^^^ 

We have introduced the Larmor frequency ujl = fisB, where fiB and m are the Bohr 
magneton and electron mass, respectively. Here, H — fiN is the pairing hamiltonian (0), 
Lz the z -component of the orbital angular momentum, the z -component of the spin, 
V the volume and S the entropy of the system. The single electron levels are Zeeman-split 
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by the magnetic field, 



et{B)=et + ujL{^t + 2T.t), (23) 



~e, = ef^ + '^{t\x^ + y^\t). (24) 

The ^-projections of the orbital and spin momentum are denoted by and Sj. The 
gyromagnetic factor of the electron is set to 2 and ef^ = 6t{B = 0). 
For spherical systems 

{t\x^ + y^\t) = '^{t\r^\t) 
and the expression for Q can be reduced to the following form: 
n = 2j2{et-fi-et) + AyG-2Tj2 log [(1 + exp(-Ei/T))(l + exp(-Ef/T))] , (25) 

t>0 t,A>0 

Et = et + uJiiAt + 2St), Et = - ujL{At + 2Si), 



et = ^{et - + A2, A = 2GY. ^Ml - ft - fi 



tA>0 



ft={l + exp{Et/T))-\ fi={l + exp{Et/T))-\ 

For normal grains {G = 0) the first two terms in Eq.(p5D vanish and Et is replaced by 
it — /i. Inserting (^) and ( P^D into ([2^), one obtains the grand canonical susceptibility of 
a spherical grain, 

X = Xd + Xp, 

as a sum of a diamagnetic and a paramagnetic contribution, xd and xp, respectively. 
The diamagnetic contribution is given by 

= E I I t)nt. (26) 

'J'^"' t,A>0 

The occupation numbers rit are for G > 

^*=2{l-^^;(or^(l-2/t)} (27) 



r) = V(4°^-//)2 + A2, /, = (l + exp{er/T})-^ (28) 
and for G = 

rit = fu ft = {l + exp{(4°) - /i)/T}"\ (29) 

The matrix element in (^) can be straightforwardly calculated with Bessel functions, 
which are eigenf unctions of our model, 

, , 5 , , A (/ + 1/2)2-1 ^ ^ 

12 



Here, k is the wave number of state t. In the single shell approximation the pairing acts 
only in the last shell. Since (t | | t) is constant within one shell, r"^ oc N and equal to its 
value without pairing. Hence, xd{G,T) is given by calculated for G = 0. 
The paramagnetic contribution, 

XP = ^ E - /*)(A* + 2S,)^ = W E - /*)(A? + 1), (31) 

is very sensitive to the pairing and the shell structure. Again, we consider small temper- 
atures (T < Tc), such that the level spacing near the Fermi level is much larger than T, 
and calculate xp by means of the single shell model. Taking into account that for the level 
with the number of magnetic substates M = 2/ + 1, 



EA 

we obtain for xp 



2 _ /(/ + 1)(2/ + 1) _ (M^ - 1)M 
~ 3 ~ 12 



XP{G = 0, T) = ^n{l - n)M{M^ + 11), n = iV,,/2M, (32) 

XP{G > 0, T) = ^n{l - n)M{M' + 11) [!"( ^ )1 (33) 

As seen, the single shell model expressions for xp can be written such that the influence 
of the pairing is expressed by a separate factor: 

Xp{G>0,T) = xpiG = 0,T)\l-(^^] ). (34) 



The temperature dependence of x is displayed in Fig. 10, which shows the case = 
3371, corresponding to the middle of the shell, where xp has its maximum. For T —>■ the 
paramagnetic susceptibility goes to zero. Hence, the grain becomes an ideal diamagnet. 
The reason is the same as for an atom. The first excited state is at the energy 2E{0) ^ 
T. i. e. only the diamagnetic part ( p6D contributes. The susceptibility for the non- 
superconducting system Xp{G = 0,T) diverges, because it costs no energy to occupy the 
magnetic substates such that the magnetic moment is finite. Hence, at small temperatures 
the susceptibility of superconducting grains is negative. It increases with the temperature, 
changes sign and reaches its maximum at T = T^. Then it decreases again proportional to 
1/T. 

At low temperatures T ~ (0.1 -=-0.2)Tc, the susceptibility for the unpaired state, xiG = 
0), takes large paramagnetic values in the open shells [^, This is an example of the 
general appearance of paramagnetism in a confined electron system, which is reviewed e. 



g. in 113. Evidence for this paramagnetic enhancement in mesoscopic normal systems 



has been foimd [29|. In our case x{G = 0) is about 10 times higher than x(G > 0). 
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Thus, measuring x these temperatures can give information whether a grain is normal 
or superconducting. 

Averaged values of x as a function of N^^'^ are shown in Fig. 11. Both cases (7 = 
and G > are displayed for T = 3. The averaged values of xiG = 0) fluctuate around 
some constant. At T = 3K, which is higher than Tc{bulk), many of grains in the range 
10^ < N < 10^ are superconductors. Hence, the averaged values of x{G > 0) are less 
than x{G = 0). For > 10^ the susceptibility x{G > 0) approaches x{G = 0). The 
finite size effects, which make the small grains superconducting at T = 3K, are no longer 
strong enough to sustain the superconducting state. For N^^^ < 15, the two different A^ 
dependences of G are reflected by the susceptibility: if G grows with decreasing A^ according 
to (0), X decreases (curve 2). If it decreases according to (0), x increases and reaches 
x{G = 0). Therefore, the measurement of the A^-dependence of the susceptibility at low 
temperatures (T < Tc) could give valuable information concerning the A^-dependence of 
G. 

Let us now derive the susceptibility for the canonical ensemble. In the presence of a 
magnetic field B the degeneracy of states with a definite seniority s is lifted, i. e. 

SE{B, s, L, S, A, S) = ujl{As + 2S,) + ^{sLS \ -r^ \ sLS). (35) 

As before, Lg, Sg, As, are the orbital and spin momenta and their projections, respec- 
tively. The subscript s indicates that these momenta correspond to a fixed seniority s and 
that they are chosen to be consistent with the Pauli principle. 

In the single shell approximation, xd is given by (^) calculated for G = 0. The reason 
is the same as for the grand canonical ensemble. The paramagnetic contribution can be 
written in a G-independent form like the canonical heat capacity if we introduce the BCS 
critical temperature Tc as the unit of the temperature and energy, 

XP = ^T.-wh'-^)Rs, (36) 
Rs= (As + 2S,)2, (37) 

= Ei^'^/T,, t = T/T,, (38) 

To calculate the canonical xp oiis needs the values of the orbital and spin momenta at a 
given seniority. This problem is solved in the Appendix. 

The results for the shells with / = 2 and / = 5 are shown in Fig. 12 in comparison with the 
grand canonical calculations. The relationship (0), which is exact for the grand canonical 
ensemble, holds with high accuracy also in the case of the canonical ensemble if the grand 
canonical A is replaced by the canonical Acan- It permits a transparent interpretation of 
the modification of the susceptibility due to the conservation of the particle number. The 
consequences of the superconductivity are expressed by the factor 1— (Acan(T)/A(0))^. The 
canonical approach washes out the sharp boundary between superconducting and normal 
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states. As illustrated in fig. 7, the canonical gap disappears gradually. Correspondingly, 
the transition from superconducting and normal values of the susceptibility is smoothed. 
The exact calculation supports this conclusion. 



7 Conclusions 

We have studied the dependence of the pair correlations on the particle number in Al 
nanometer-scale grains, 10^ < N < 10^. Using the spherically symmetric infinite well as a 
model of the field, we have calculated the gap parameter A as functions of A^. Comparing 
it with experimental data from tunneling experiments, the pairing coupling constant G{N) 
is fixed as function of the number of electrons A^. 

In our model, the bunching of the electron levels (shell structure) strongly enhances the 
pair correlations in the small grains. This enhancement is so strong that fitted coupling 
constant A(A^) = G{N)pf decreases with decreasing A^, in order to account for the more 
modest increase seen in experiment. We interpret this as a consequence that our spherical 
model has too pronounced a shell structure. The bunching of the electronic levels in realistic 
grains is most likely weaker due to deviations of the shape from the ideal sphere, surface 
roughness and impurities. Thus, the enhancement of the superconductivity due to the shell 
structure is expected to be weaker and a different function A(A^) will fit the experimental 
data on A(A^), which for the limiting case of no shell structure increases with decreasing 
N. 

The enhancement of the pair correlations by the shell structure should be considered as 
a mechanism that exists in addition to the increase of the effective pair coupling constant A 
caused by the modification of the phonon spectrum in small grains. The completely differ- 
ent A^ dependence of A found for the spherical and structureless models demonstrates that 
a careful estimate of the shell structure in realistic grains is needed in order to determine 
the A^ - dependence coupling constant from the data. 

The averaged magnetic susceptibility for the spherical grains in the normal state turns 
out to be strongly paramagnetic. Evidence for this paramagnetic enhancement in meso- 



scopic normal systems has been found [2£]. In our spherical model the superconductivity 
strongly reduces the paramagnetism, but the susceptibihty remains large and positive. 
Thus, the surprising prediction is that small superconductors may be paramagnetic. In 
realistic grains with a less pronounced shell structure the pair correlations may be strong 
enough to make the susceptibility negative. Our study shows that the susceptibility of 
grains composed of superconducting material results from an intricate interplay between 
the pair correlations and the spatial confinement of the electrons. 

The properties of the superconducting grains with the number of atoms below 10^ 
are significantly modified by the fixed number of electrons. Instead of a sharp phase 
transition an extended transition region in temperature appears. A transition temperature 
can be defined as the value where the the specific heat shows the most rapid drop. This 
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critical temperature agrees rather well with Tc calculated in the grand canonical BCS 
approximation and can be adopted as the transition temperature Tc for small systems. 

Part of this work was performed during the visit of two authors (N.K. and V.M.) to 
Forschungszentrum Rossendorf. They thank Prof. Prade for hospitality. This work is 
supported by the INTAS (grant INTAS-93-151-EXT). 
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8 Appendix 



The sum Rg appearing in Eq.(|37|) amounts to the difference of the sums of squares of 
orbital and spin momentum projections corresponding to states with N = s and = s — 2 
particles, 

Rs = Rn=s — Rn=s-2 (39) 



RN = T.{i^niL.+2s.r\i^n 

i 

Rjy includes the expectation values of the z-projection orbital and spin momentum opera- 
tors (1/2 and Sz). The total number of states ipf at a given N is 

M = 21 + 1, 

with I being the orbital momentum. In what follows we assume that is even and 
equal to 2k. As known the exclusion principle requires that orbital momentum and 




spin functions entering into ipi are basic vectors of conjugate or dual representations of the 
permutation group and their Young diagrams should correspond to each other by exchange 
of rows and columns. The diagrams of spin functions consist of k + S squares in the first 
row and k — S in the second one. S is the spin taking the only value for each diagram and 
ranging from to k = N/2 for a fixed even integer A^. The orbital permutation symmetry 
is characterized by two column diagrams, the whose length of columns is equal to A; + S" 
and k — S. In general, such diagrams involve several states with orbital momenta L^g. We 
assume that ipi^ is a direct product of an orbital function ip{k, S, L^s, and a spin one 
x{S, S) where A and S are eigenvalues of the operators and 5*2 respectively. 

RN=2k= E Mk,s,Lks)x{S,^)\{Lz + 2SzY\y^{k,s,Lks)x{S,^)) (40) 

As the sums over A and S are independent, Rn can be divided into orbital and spin parts: 

Rn = Rn + (41) 
<=2fc = E s, Lks) I Ll I ^{k, s, Lks)){2S + 1), (42) 

LkS,S 

A k 

Rl=2k = E(X(5,S) I ASl I x(^,S))45 = E S{S+l){2S+l)d,s, (43) 

5,2 5=1 

where dks is the dimension of the orbital space with the permutation symmetry described 
by the two column Young diagrams mentioned above. 

( M + l\ ( M + l\ 2S + 1 
^^^= k-S k + S 2MTT ^''^ 
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j:{2S + 1)45 =[2k ) 

Eq.(P5D indicates that can be computed straightforwardly whereas finding R^ needs 
the determination of the set of L^s- 

This task can be removed by taking into account that external product of two com- 
pletely antisymmetric functions ip{k + S,L,A) {k + S squares in the only column of the 
Young diagram) and (f{k — S, L', A') {k — S squares) gives rise to a series of basic vectors 
of irreducible representations with two column diagrams. Each vector arises only once and 
the lengths of columns vary from {k + S, k — S) up to {2k, 0), i.e. the decomposition of this 
external product contains ip{k, S' , L^s, A) functions with S'^^^ = S and S!^^^ = k. This 
decomposition permits to reduce Eq.(^2|) to the following 

R^=2k = i:T.{(^ik + S,L,A)\Ll\^{k + S,L,A)){ ^ )+ 

S=OL,A 

{^{k - S, L, A) I Ll I ^{k - S, L, A)) ( ^ ) } (2 - ^5,0) (46) 

The orbital momenta L oiip have to be compatible with antisymmetry of these states. This 
points out an elementary way to compute them by summing single particle projections 
among which should not be identical ones. 
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Figure 1: Averaged values of AbuikiT = 0)/Ae v.s. N^^^. The pairing gap of bulk Al is 
AbuikiT = 0) and Ae = ep+i — ep. ep, ep+i are energies of the Fermi level and the next 
higher one. The meaning of the averaging is explained at the end of section ^ 
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Figure 2: Pairing gap A{T,N) v.s. N^^^. Dashed lines: T = 0, solid lines: T = 3 K. 
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Figure 3: Tc/E{0) v.s. the occupation degree of the shell n = Nsh/2M. 
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Figure 4: Averaged values of 2E{0)/Tc and 2A(0)/Tc v.s. N^/^. 
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Figure 5: (a) v.s. N^^^. Points represent experimental data |^ on A, which are 
converted by mean of eq. (|12|) into T^. The sohd hne is the averaged calculated with 
G = 1.94:N~^'^'^ meV. (b) Averaged pairing gap v.s. A^^/^ calculated at T = 0, 3 and 5 
K. The horizontal line gives A(0) of bulk Al. 
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Figure 6: (a) Tc v.s. N^^^. Points represent experimental data |^ on A, which are 
converted by mean of eq. into Tf.. The sohd hne is the averaged calculated with 
G = 3.21A^-iexp(-25A^-°-26) eV. (b) See the caption in §b. 
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Figure 7: Pairing gaps calculated v.s. T/T^ at different values of M for the half filled 
shell. Solid and dashed lines correspond to canonical and BCS results, respectively. 
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Figure 8: Specific heat C of the single shell v.s. T /Tc at different M. Solid and dashed 
lines correspond to canonical and BCS results, respectively. 
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Figure 9: Susceptibility (in units | xl \ -10^ where xl is the Landau expression for 
diamagnetism of the degenerate free electron gas) v.s. N^^^. Dashed lines: G = 0, solid 
lines: G > 0. 
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Figure 10: Susceptibility (in units \ xd \ {G 
lines give x at G = and G > 0, respectively. 



0,T = 0)) v.s. T/Tc. Dashed and solid 
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Figure 11: Averaged susceptibility (in units | xl \ -10^) v.s. N^^'^ at T = 3 K. Curves 
1,2,3 correspond to G = 0; G = 1.94iV-o-^^ meV] G = 3.21A^-i exp(-25A^-°'26) eV, 
respectively. 
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Figure 12: Grand canonical and canonical paramagnetic susceptibilities (in units Xp{Cr = 
0,T = Tc) ) v.s. T/Tc. Curves 1,2 represent the results of canonical calculations of 
XpiG > 0) for the shell with M = 11, N^h = 10 (curve 1) and for the shell with M = 5, 
Nsh = 5 (curve 2). 
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